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ABSTRACT 



Context. Observations indicate that the 'quiet' solar photosphere outside active regions contains considerable amounts 
of magnetic energy and magnetic flux, with mixed polarity on small scales. The origin of this flux is unclear. 
Aims. We test whether local dynamo action of the near-surface convection (granulation) can generate a significant 
contribution to the observed magnetic fiux. 

Methods. We have carried out MHD simulations of solar surface convection, including the effects of strong stratification, 
compressibility, partial ionization, radiative transfer, as well as an open lower boundary. 

Results. Exponential growth of a weak magnetic seed field (with vanishing net flux through the computational box) is 
found in a simulation run with a magnetic Reynolds number of about 2600. The magnetic energy approaches saturation 
at a level of a few percent of the total kinetic energy of the convective motions. Near the visible solar surface, the 
(unsigned) magnetic flux density reaches at least a value of about 25 G. 

Conclusions. A realistic flow topology of stratified, compressible, non-helical surface convection without enforced recir- 
culation is capable of turbulent local dynamo action near the solar surface. 
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1. Introduction 

Magnetic fields in the quiet solar photosphere (often 
referred to as 'internetwork fields') are of considerable 
interest in connection with the heating of the upper solar 
atmosphere and as a possible example of local fast dynamo 
action. Owing to their apparently 'turbulent' nature, 
with mixed polarity on small spatial scales, their extended 
range of field strengths, and their weak polarization signals, 
quiet-Sun (QS) magnetic fields are a difficult observational 
target. Quantitative results from observations based upon 
different methods do not yet provide a fully consistent pic- 
ture ( e.g., [Lin & Rimmelc, 1999; Sanchez Almeida & Lite^, 
,2000; K homenko et al.l.l2003tlSanchez Almeida et al.l.l2003|; 
Domfnguez Cerdena et al. ~ l2003l: iLites fc Socas- Navarre 



2001" 



iDominguez Cerdena et al.[ 



I2006all]: 



Truiillo Bueno et al.l . 12004 ISanchez Almeidal . l2003l 



20051) ■ but it appears to be established that small-scale, 



mixed-polarity magnetic fields are ubiquitous in the quiet 
Sun and contribute significantly to the total magnetic 
energy and unsigned flux in the photosphere outside active 
regions. 

The origin of these small-scale fields is not yet fully 
clarified. It seems plausible that rise and emergence of 
magnetic flux from the deep convection zone in the form 
of small bipoles as well as the debris from decaying ac- 
tive regions contribute to the QS magnetic fields. On the 
other hand, th e increasing flux replen ishment rates towards 
smaller scales ([Hagenaar et al.l I2003D and the absence of a 
significant variation of the QS fiux between solar minimum 
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and maximum (iTrujillo Bueno et al.l . l2004[ ) indicate that a 
considerable part of the QS magnetic fields may be pro- 
duced by a self-excited turbulent dyn amo operating in the 
near- surface layers of the Sun (e.g., iPetrovav fc Szakalvl 
Il993l) . Since the turnover time of the corresponding dom- 
inant convection pattern (granulation) of about 10 min is 
much smaller than the solar rotation period, this flow is 
practically unaffected by the Coriolis force and thus ex- 
hibits no net helicity. The principle of such a non-helical 
surface dynamo has been demonstrated in incompressible 
closed-box MHD s imulations of thermal convection at hig h 
Rayleigh number (jCattaneol . 119991 : ICattaneo et al.l . |2003[ ). 
The question is whether these results carry over to the 
case of realistic solar granulation. This is not at all ob- 
vious since numerical simulations of solar surface convec- 
tion show that the powerful surface cooling by radiation 
and the strong stratification prevent a significant recircu- 
lation of the downf l qwing plasma in the near-surface layers 
(iStein fc Nordlundl . (19891) . As a consequence, small-scale 
magnetic flux is rapidly pumped into the deeper layers of 
the convection zone and is th us lost for further amplifi ca- 
tion by the near- surface flows (jStein fc Nordlundl . [20031) . In 
a simulation with artificially closed boundaries, such pump- 
ing is suppressed and a strong local recirculation main- 
tained, so that a local dynamo is efficiently sustained. 



In this paper, we present the first example of dynamo 
action in a realistic simulation of solar surface convection 
with an open lower boundary. Our results indicate that a 
realistic fiow topology of strongly stratified, compressible, 
and non-helical surface convection without enforced recir- 
culation is capable of driving a turbulent local dynamo. 



2 



A. Vogler and M. Schiissler: A solar surface dynamo 



Table 1. Simulation parameters. The estimate of the mag- 
netic Reynolds number is based on the rms flow velocity 
and a length scale of 1 Mm (granulation scale) . Az and Ax 
are the vertical and horizontal grid spacing, respectively. 



Run 


height /width 


Az/Ax 


V 


Rm 




[Mm] 


[km] 


[lO^Ocm^s-i] 


(approx.) 


A 


1.4/6.0 


14/20.8 


11.1 


300 


B 


1.4/6.0 


14/10.4 


2.5 


1300 


C 


1.4/4.86 


10/7.5 


1.25 


2600 



2. Numerical model 



We have used the MURaM code (IVogler et all . l2005t 
iVoded . [200l to carry out local-box MHD simulations of 
solar surface convection with grey radiative transfer. The 
computational domain covers the height range between 
about 800 km below and 600 km above the average height 
of the visible solar surface and has a horizontal extension 
between about 5 Mm and 6 Mm (see Tab.[T|). 

The side boundaries are periodic in both horizontal di- 
rections. The boundary condition at the bottom {z = 0) 
permits free in- and outflow of matter: the upflows are as- 
sumed to be vertical, = Vy = , dzVz — 0, in the down- 
flows the vertical gradients are set to zero, dzVx — dzVy — 
dzVz = 0. The upper boundary is closed for the flow. The 
magnetic field is assumed to be vertical at the upper and 
lower boundaries: = By = , dzB^ = 0. F urther details 
of the boundary conditions are described in IVdgler et al.l 
(|2005() . Horizontal fields which get carried towards the lower 
boundary in downflows nevertheless leave the simulation 
domain by means of magnetic diffusion across the bound- 
ary: at any thime, the diffusive boundary layer at the bot- 
tom adjusts its width such that the resulting diffusive flux 
matches the incoming advection flux into the boundary 
layer. Since the magnetic field strength in the lower part 
of the box always remains significantly below the equipar- 
tition value with the convective flows, it is thus guaranteed 
that the magnetic flux carried by the downflows leaves the 
box unimpeded. In order to always ensure that the diffusive 
boundary layer is well resolved numerically and to reduce 
the effect of the numerical boundary on the magnetic fields 
in the bulk of the simulation domain, the magnetic diffusiv- 
ity is increased in a region of 150 km thickness at the bot- 
tom boundary. The lower boundary condition also prevents 
the advection of horizontal magnetic flux into the box from 
below. In the real Sun, of course, magnetic flux from the 
deeper layers is probably advected into the surface layers 
and may influence the properties of the QS magnetic field. 
In our experiment, we intentionally exclude this source of 
magnetic fiux in order to study local dynamo action in iso- 
lation. 

While we have a constant value of the magnetic dif- 
fusivity in the box (outside the region of enhanced diffu- 
sivity near the bottom), the code uses an artificial viscos- 
ity that varies in space and time (see Voglcr et al., 2005). 
Consequently, we can evaluate the magnetic Reynolds num- 
ber, for a simulation run, but a quantitative esti- 
mate of the hydrodynamic Reynolds number. Re, and of 
the magnetic Prandtl number, Pr^ = Rax/ Re, is diffi- 
cult to obtain. Since both, the explicit magnetic diffusiv- 
ity and the artificial viscosity, lead to a diffusive cutoff 
at the scale of the grid resolution. Re and i?ni are of the 




time 



Fig. 1. Total magnetic energy (normalized to the energy 
of the initial seed field at f = 0) in the simulation box 
as a function of simulated solar time for the three runs 
specified in Tab. [TJ All runs show an initial rapid energy 
increase due to flux expulsion of the seed field by granu- 
lation. Thereafter, run A shows an exponential decay with 
an e-folding time of roughly one hour, while run B is ap- 
proximately marginal at a low energy level. Run C exhibits 
exponential growth with a time scale of about 10 minutes 
and approaches a saturation level of a few percent of the 
total kinetic energy of the convective flow. 



same order of magnitude and thus is about unity at 
the smallest resolved spatial scales. For smooth flows on 
larger scales, the artificial viscosity is much smaller than 
the magnetic diffusivity, so that our 'effective' magnetic 
Prandtl number on these scales is smaller than unity. Given 
these conditions, our results do not provide a proper ba- 
sis for commenting on the much debated question of the 
dependenc e of turbulent dynarnq action on the value of 
Prm ('e.g.. iBoldvrev k. Cattaneol. 12004 iPontv et all 12005 : 



Sch ekochihiin et al.1 120051 : iBrandenburg fc Subramanian . 
,20051) . 

We have carried out three simulation runs with different 
magnetic Reynolds numbers, For all runs, a weak mag- 
netic seed field with zero net fiux was introduced into a fully 
developed nonmagnetic convection pattern. The seed field 
was purely vertical, with polarity variations corresponding 
to a checkerboard-like 4x4 horizontal planform and a con- 
stant field strength of |i3o| = 10 mG. The box dimensions, 
grid resolution, and magnetic diffusivities for the runs are 
given in Tab. [TJ 

3. Results 

Figure [1] shows the magnetic energy for the three simu- 
lation runs as a function of time. In run A {Rm ~ 300), 
the magnetic energy decays exponentially after a brief ini- 
tial amplification due to flux expulsion acting on the initial 
magnetic configuration. Run B {R^n ~ 1300) appears to be 
close to the point of marginal dynamo excitation. The mag- 
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Fig. 3. Probability density function (PDF) for the vertical 
field component at three different geometrical height levels. 
Dashed curve: z = —370 km (about 450 km below the visi- 
ble surface); solid curve: z = 80 km (roughly corresponding 
to the average level of tr — 1); dotted curve: z = 400 km 
(about the average level of tr = 0.01). Shown are time 
averages over about 20 minutes around t ~ 4.5 h. 



netic energy levels out well within the kinematic regime, 
with maximum field strengths several orders of magnitude 
below local equipartition with respect to the kinetic energy 
density of the flow. Finally, run C (i?m ~ 2600) exhibits 
exponential growth with an e-folding timescale of roughly 
10 minutes. After about three hours of simulated solar time, 
the magnetic energy approaches a saturated level of about 
3% of the total kinetic energy of the convective flow. 

For a snapshot during the saturation phase of run C, 
Fig. [5] shows maps of the (bolometric) brightness and of 
the vertical magnetic field on two surfaces of constant 
(Rosseland) optical depth, tr. Around tr = 1, the field 
exhibits an intricate small-scale mixed-polarity structure, 
which extends down to the diffusive length scale. The (un- 
signed) mean magnetic fiux density at this level has reached 
a value of about 25 G. At the surface tr = 0.01, about 
300 km higher in the atmosphere, the spatial distribution 
of the field is significantly smoother and more closely asso- 
ciated to the intergranular downflow lanes. The mean flux 
density has decreased to about 3 G. This strong decrease 
indicates the absence of significant dynamo driving in the 
convectively stable layers above tr — 1, so that the field 
decays rapidly with height, owing to its small horizontal 
spatial scale near the visible solar surface. 

Figure [3] shows the average probability density function 
(PDF) of the vertical magnetic field, determined during the 
saturation phase of the dynamo, at three height levels. The 
PDFs have the form of stretched exponentials, indicating 
a strong intermittency of the magnetic field at all heights. 
The strongest magnetic features occasionally reach vertical 
field strengths beyond 1 kG near tr = 1. 

Energy spectra for the vertical components of the near- 
surface magnetic field and velocity as a function of hori- 
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Fig. 4. Energy spectra based on the vertical components 
of velocity and magnetic field, respectively, as functions of 
horizontal wave number, k\^. The values are taken at 2; = 0, 
corresponding to a depth of about 80 km below the aver- 
age average level of tr = 1), at i ~ 5 h. The kinetic energy 
spectrum (solid curve) peaks at ~ 3...4Mm~^, roughly 
corresponding to the typical scale of granules. The mag- 
netic energy spectrum has a broad maximum around wave 
numbers of about 30Mm~^, corresponding to length scales 
of at least an order of magnitude smaller. 



zontal wave number, /ch, are given in Fig. [51 The spectral 
magnetic energy shows a broad peak at fch — 30, which 
corresponds to a wavelength of about 200 km. At the high- 
wavenumber end of the spectra, the magnetic and kinetic 
energies become less disparate. The remaining deviation 
from equipartition is due to the anisotropy resulting from 
the strong stratification. 

We find that convective downward pumping of fiux in 
fact has a significant effe ct on the energy ba l ance o f the dy- 
namo, as conjectured bv I Stein fc Nordlundl (|2003l ). At any 
given height, the time dependence of the horizontally av- 
eraged magnetic energy density, emag, is governed by the 
equation 9femag = Wl - W^j - dV^/dz. W^L is the rate 
of work against the Lorentz force, Wj is the Joule heating 
rate, and Vz is the vertical component of the Poynting flux, 
the advective part of which measures the draining of mag- 
netic energy due to convective pumping. All quantities are 
meant to be horizontal averages. The advective Poynting 
fiux is negative throughout the convectively unstable parts 
of the simulation domain, confirming that any growth of 
magnetic energy in the system must have its source inside 
the domain. The diffusive part of Vz is found to be negli- 
gible in the convecting layer. In the absence of convective 
pumping, the difference Wi^ — Wj would be a measure for 
the growth of the magnetic energy during the exponential 
growth phase. In our case, more than 80% of this difference 
is indeed carried downwards by means of the term —dVz/dz 
and leaves the box through the bottom boundary. However, 
the effect only reduces the growth rate but does not shut 
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Fig. 2. Snapshot from the dynamo rnn C, taken about 5 honrs after introducing the seed field. The vertically emerging 
bolometric intensity (brightness, left panel) reveals a normal solar granulation pattern. The other panels show the vertical 
component of the magnetic field on two surfaces of constant (Rosseland) optical depth, tr. Near the visible surface (middle 
panel, tr, = 1, grey scale saturating at ±250 G), the magnetic field shows an intricate small-scale pattern with rapid 
polarity changes and an unsigned average flux density of 25.1 G. About 300 km higher, at the surface tr = 0.01 (right 
panel, grey scale saturating at ±50 G), the unsigned average flux density has decreased to 3.2 G and the field distribution 
has become considerably smoother, roughly outlining the network of intergranular downflow lanes (darker areas on the 
left panel). 



down dynamo action if the magnetic Reynolds number is 
sufficiently large. 

4. Discussion 

Our main finding is that a realistic flow topology of strongly 
stratified convection in the near-surface layers of the Sun 
is capable of sustaining dynamo action. Downward pump- 
ing in an open box has a signiflcant impact on the energy 
balance, but is not able to shut down the dynamo, presum- 
ably because there is sufficient local recirculation (e.g. by 
turbulent entrainment of downflowing material into upflow 
regions) to amplify the magnetic field near the surface. A 
more detailed analysis will have to be carried out in order 
to clarify the physical mechanism at work here. Strictly 
speaking, a dynamo simulation would have to be run for 
several global diffusive timescales before transients can be 
definitively ruled out. This is practically not feasible in our 
case since the global diffusion time for run C is about 18 
days. On the other hand, the exponential amplification of a 
miniscule seed field to substantial saturation levels of flux 
over several hours, with a growth time corresponding to 
the granulation time scale, lends credibility to physical rel- 
evance to the resiilts. 

While demonstrating the possibility of local dynamo ac- 
tion by granulation in principle, the quantitative results of 
our simulations should be interpreted with caution. Owing 
to our subgrid model for the viscosity, the magnetic Prandtl 
number cannot be uniquely defined. We have argued that 
its value is significantly smaller than unity at the spatial 
scales where the magnetic energy peaks, but certainly our 
numerical experiment is not suited to address the question 
of turbulent dynamo action in the limit of very small Pri„. 
Likewise, the saturation level of our dynamo run proba- 
bly depends on the magnetic Reynolds number (as well as 
on Prm), so that a quantitative comparison with observa- 
tional results in terms of mean flux densities and PDFs 
would certainly be premature. Such comparison would also 



have to take into account, in addition to the self-excited 
dynamo, the advection of flux from below and the process- 
ing of magnetic debris from decaying active regions. Still, 
we deem it noteworthy that levels of magnetic flux of the 
observed order of magnitude can, in principle, be produced 
by simulations of the kind presented here. 

This is just a first step and much remains to be done 
in order to establish quantitatively reliable results of near- 
surface local dynamo action. This includes studying the ef- 
fects of deeper and wider computational boxes as well vari- 
ations of the boundary conditions. Clearly, the dependence 
of the saturation level on the magnetic Reynolds and mag- 
netic Prandtl numbers has to be investigated through series 
of controlled numerical experiments. Comparison with ob- 
servational results will be crucial in order to evaluate the 
contributions by flux advection from the deeper layers. 

Acknowledgements. We are grateful to Robert Cameron and 
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